Simulations of a classical spin system with competing superexchange and 

double-exchange interactions 
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Monte-Carlo simulations and ground-state calculations have been used to map out the phase di- 
agram of a system of classical spins, on a simple cubic lattice, where nearest-neighbor pairs of spins 
are coupled via competing antiferromagnetic superexchange and ferromagnetic double-exchange in- 
teractions. For a certain range of parameters, this model is relevant for some magnetic materials, 
such as doped manganites, which exhibit the remarkable colossal magnetoresistance effect. The 
phase diagram includes two regions in which the two sublattice magnetizations differ in magnitude. 
Spin-dynamics simulations have been used to compute the time- and space-displaced spin-spin cor- 
relation functions, and their Fourier transforms, which yield the dynamic structure factor S{q, uj) 
for this system. Effects of the double-exchange interaction on the dispersion curves are shown. 
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I. INTRODUCTION 



Doped rare-earth manganites with general chemical 
formula Rei-^^AajMnOa (where Re denotes a rare-earth 
and A is typically Ca or Sr) have been shown to exhibit 
colossal magnetoresistance, and a wide range of other 
physical properties, which are strongly dependent upon 
doping yj. The physical mechanisms responsible for the 
unusual magnetoresistance properties of these materials 
have not been understood yet, and it has been suggested 
that double-exchange interactions between Mn+^ and 
Mn+^ ions are important. However, there is much ongo- 
ing debate as to whether double-exchange alone can in- 
deed give rise to colossal magnetoresistance or if further 
mechanisms have to be considered [J. It is thus impor- 
tant to further understand properties associated with the 
double-exchange theory. 

Equal valence manganese ions in the parent com- 
pounds interact via a superexchange mechanism, whereas 
the different valence manganese ions, introduced through 
doping, are coupled via a double-exchange interaction 
QiM @| ■ We investigate a simple model of classical spins 
in which nearest-neighbor spins interact via a compet- 
ing antiferromagnetic superexchange and a ferromagnetic 
double-exchange interaction, where the ratio of these two 
interactions mimic the effects of the doping x. Proper- 
ties of the pure Heisenberg model (only superexchange 
interaction) are well-known and recently a high- 
accuracy Monte-Carlo study obtained the paramagnet- 
ferromagnet transition temperature and associated static 
critical exponents of the pure double-exchange model p| . 
Recent studies of the phase diagram of double-exchange 
systems using a different approach from ours are given 
in Ref.Q. We use a combination of Monte-Carlo sim- 
ulations and ground-state calculations to determine the 
phase diagram of the model with varying ratios of the su- 
perexchange and double-exchange interaction strengths. 
We also use spin-dynamics techniques to study the dy- 
namic structure factor and obtain the dispersion curve of 
the double-exchange model. 



II. MODEL AND METHODS 

The model considered here can be described by the 
Hamiltonian 

H — ~JaF ^ Sr-Sr'— J_D_E ^ v^l + Sr • Sr' , (1) 
<r,r'> <r,r'> 

where Sr = {Sr^ , Sr^ , S/') is a three-dimensional clas- 
sical spin of unit length at site r. The first term, with 
Jaf < 0, is the antiferromagnetic superexchange inter- 
action and the second term describes a ferromagnetic 
double-exchange term with coupling constant Joe > ^ 
between nearest-neighbor pairs of spins. We consider 
L X L X L simple cubic lattices with periodic bound- 
ary conditions, and divide the system into two sublat- 
tices, A and B, with magnetizations and ms, re- 
spectively. The scalar product is defined as niA • m_B = 
|m^||mB|cos0 = y, where clearly — 1 < y < 1. 

In the ground state, the condition of energy minimiza- 
tion leads to the expression 



y^iJDE)ViMJAl 



1, 



(2) 



with \Jaf\/Jde > l/\/8 = 0.35355..., and for 
\Jaf\/Jde < 0.35355... the ground state is ferromag- 
netic. Eq.© admits the possibility of non-colinear 
phases and/or sublattice magnetizations with different 
magnitudes. Note that y = for \Jaf\/ Jde = 1/2, 
which means that 9 — tt/2 and/or that the magnetiza- 
tion in at least one sublattice is zero in this case. 

For T > we use Monte-Carlo simulations to mea- 
sure thermodynamic quantities of the model in order to 
determine the phase diagram. We used the Metropo- 
lis algorithm, lattice sizes L = 12 and 24, and typically 
we discarded 10000 Monte-Carlo steps for thermalization 
and used 10000-30000 steps in computing averages. 

The dynamics of the spins are governed by the equa- 
tions of motion 



dt 



Sr = Sr X (-Vr^), 



(3) 



2 



where the effective field 



1.0 



Jdi 



E 



x/1 + Sr • Sr' 



(4) 



and the sums are over r' nearest-neighbor to r. These 
equations were solved numerically using a method based 
on second-order Suzuki- Trotter decomposition of expo- 
nential operators to a maximum time of tmax = 
440/Jb£. 

The dynamic structure factor S'(q, w) for momentum 
transfer q and frequency transfer w, observable in neu- 
tron scattering experiments, is given by 

/ + 00 I. 
expM)C'=(r-r',t)^, 
rr- V^TT 

(5) 

where C'^(r — r',t) is the space-displaced, time-displaced 
spin-spin correlation function defined, with k — x,y, or z, 
as C'=(r-r',t) = (5r'^'(05r''^'(0)) - (5'r'' W)(^r'^'(0)). The 
displacement r is in units of the lattice unit cell length 
a. More details on the spin-dynamics methods are given 
inRefs. U^E^. 



III. RESULTS 

Our preliminary results show that \m.A\ — \m.B\ at 
all temperatures for \Jaf\/ Jde < I/a/S, decreasing 
smoothly from |m^| = Imsl = 1 at low T, towards 
zero as T increases. For \Jaf\/Jde > l/v^ we find 
that at low temperatures |m^| ^ Imsj, as illustrated 
in Fig. 1 for \Jaf\/Jde = 0.9. For this range of 
interaction strengths our results indicate that \mA\ = 
\{JDE?/i'i\JAF\^) - 1| and Imsl = 1 at T = 0, and 
Eq.(j2l then implies that = for \Jaf\IJde < 1/2, 
9 = TT for \Jaf\/Jde > 1/2, and the sublattice mag- 
netizations reduce to |m^| = and Imsl — 1 for 
\Jaf\/ Jde — 1/2. The phase-transition temperatures 
are obtained from analyses of the sublattice magnetiza- 
tions for L = 24, hence there are some uncertainties in 
these estimates due to finite-size effects. 

The phase diagram that we obtained from the anal- 
ysis of the sublattice magnetizations, shown in Fig. 2, 
has five different regions: (i) ferromagnetic (F) for low 
T and small values of \ Jaf\/ Jde] (ii) antiferromagnetic 
(AF) for moderate T and high values of \ Jaf\/ Jde', (hi) 
paramagnetic (PM) at high T; and two low tempera- 
ture phases, which we label as (iv) region I for l/y/S < 
\Jaf\/Jde < 1/2, and (v) region II for \Jaf\/Jde > 
1/2. Regions (i)-(iii) are characterized by \mA\ — Imsj, 
with = ms in the ferromagnetic phase and niA — 
— rns in the antiferromagnetic phase. In contrast, in re- 
gions I and II we find that Im^l ^ Ims | and although our 
simulations suggest that 9 = and tt in regions I and II, 
respectively, we cannot rule out possible small deviations 
from the aligned and anti-aligned sublattice spin config- 
urations in the respective regions. The question marks 
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FIG. 1: Sublattice magnetization for \ Jaf\/ Jde 
lattice size L = 24. 
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FIG. 2: Phase diagram. When not shown, the error bars are 
smaller than or of the size of the points and the solid lines 
simply join the data, guiding the eyes. 



(?) in the phase diagram (see Fig. 2) indicate that the 
points at which the F-PM and AF-PM transition lines 
join other phase boundaries have not been determined 
accurately yet. For the pure double-exchange ferromag- 
netic model (Jaf = 0), the critical temperature (F-PM 
transition point) is Tc^^ = 0.74515(7) Jdb 6], whereas 
for the pure Heisenberg model T^'^ = 1.442929(77)| Jaf| 
5] for either ferromagnetic or antiferromagnetic interac- 
tion. 

Fig. 3 shows a comparison of the dispersion relations 
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FIG. 3: Dispersion relations for the Hoisenberg fcrromagnet 
(pure superexchange model, labeled as SE) and the pure fer- 
romagnetic double-exchange model (DE). 

for the pure Heisenberg ferromagnet and for the pure 
double-exchange ferromagnetic model. The results at 
T = correspond to the linear spin-wave theory, and 
data at the respective critical temperatures are from our 



simulations. 



IV. CONCLUSIONS 

We studied a system of classical spins coupled via com- 
peting antifcrromagnctic superexchange and ferromag- 
netic double-exchange interactions. As expected, when 
the former is much weaker than the latter, the system 
is ferromagnetic and in the reverse case the system is 
antiferromagnetic. In addition, wc found two other low 
temperature phases, in which the sublatticc magnetiza- 
tions differ in magnitude and are close to being ahgned 
or anti-aligned, depending on the ratio of the two in- 
teraction strenghts. At high temperature the system is 
paramagnetic. Comparing the dispersion curves we see 
that spin- waves in the pure double-exchange model occur 
at lower frequency transfer than in the pure Heisenberg 
ferromagnet. 
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